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We consider one-dimensional (ID) interacting electrons beyond the Dzyaloshinskii-Larkin the- 
orem, i.e., keeping forward scattering interactions among the electrons but adding a non-linear 
correction to the electron dispersion relation. The latter generates multi-loop corrections to the 
polarization operator and electron self-energy thereby providing a variety of inelastic processes af- 
fecting equilibrium as well as non-equilibrium properties of the ID system. We first review the 
computation of equilibrium properties, e.g., the high frequency part of the dynamical structure fac- 
tor and corrections to the electron-electron scattering rate. On this basis, microscopic equilibration 
processes are identified and a qualitative estimate of the relaxation rate of thermal carriers is given. 



I. INTRODUCTION 

In systems of interacting electrons in one dimension 
certain physical phenomena such as damping of the el- 
ementary excitations and energy relaxation, are beyond 
the reach of the conventional Tomonaga-Luttinger (TL) 
model; see, e.g., Ref. [1 for a review. This may be 
understood from the fact that, for such a model, the 
spectrum of the interacting electrons is in exact one- 
to-one correspondence with the one of free bosons P In 
diagrammatic terms this correspondence is generally re- 
ferred to in the literature as the Dzyaloshinskii-Larkin 
(DL) theorem™ see also Refs. [5-8 for reviews. In- 
elastic processes responsible for a finite life-time of the 
elementary excitations and relaxation in one dimension 
may, however, be accessed by abandoning Tomonaga's 
approximation of linearizing the electron dispersion re- 
lation^ at the expense of losing the exact solvability of 
the model. A number of methods have been developed, 
especially during the last decade, to address these issues 
thereby generating a considerable amount of work; see 
Ref. [1 and references therein. The purpose of this pa- 
per is to emphasize the importance of yet another ap- 
proach based on the direct analysis of the lowest order 
multi-loop diagrams which are beyond the DL theorem. 
A review of the method for the computation of equilib- 
rium properties will first be presented. On this basis, 
microscopic processes responsible for the equilibration of 
one-dimensional (ID) electrons will be identified. The 
analysis will be qualitative but will nevertheless allow us 
to give an estimate of the relaxation rate of thermal car- 
riers in quantum wires. This estimate agrees with the 
recent literature on the subject 

II. DISPERSION NON-LINEARITY AND 
INTERACTIONS 

The TL model is exactly solvable thanks to two cru- 
cial approximations related to the shape of the elec- 
tron dispersion relation and to the nature of the scat- 
tering processes between the electrons. The former con- 
sists in linearizing the electron dispersion relation: = 



±vp(k =F kp), where ±kp are the Fermi points, ± cor- 
respond to the chiralities of the fermions (+ for right 
movers near +kp and — for left movers near —hp), and 
vp is the Fermi velocity. As will be seen shortly, it is this 
approximation which is responsible for the fact that low- 
energy excitations are free bosons, i.e., quantized sound 
waves^of the ID system. The properties of these collec- 
tive modes are encoded in the density-density correlation 
function, U(x, r) = (p(x, r)p(0, 0)), where p is the density 
operator. For non-interacting electrons, this correlation 
function reduces to a 1-loop diagram, see Fig. |2|a). In 
imaginary time and in the linear spectrum approximation 
it readP^ 

i4°W) = ±^-^ — , a) 

Z7T iuj =p Vpq 

where uj is the energy of the boson of momentum q and 
± refer to the chirality of the fermions. The total po- 
larization operator is given by: n(°) =nf +nL° } . The 
imaginary part of the retarded polarization operator de- 
fines the dynamical structure factor (DSF): S(uj,q) = 
— ^ ^sIl R (uj , q) . Combining this definition with Eq. 
yields 

S^(u,q) = ^S(w-v F \q\), (2) 

7T 

where we assume that uj > for simplicity. Equation ([2| 
together with the /-sum rule^ 

roc q 2 

/ duj uj S(uj,q) = n , (3) 

Jo 2m 

where n is the average electron density along the wire, 
show that indeed the free boson of energy uj = vp\q\ 
exhausts the /-sum rule. 

A coupling between the collective modes can be 
achieved by taking into account small deviations from 
linear dispersion. Away from half-filling, the first cor- 
rection is due to the curvature of the electron dispersion 
relation: = ±vp(k=fkp) + (k^fkp) 2 /2m, which intro- 
duces an additional parameter, m, the band mass such 
that: kp « mvp, at small fillings. Including this non- 
linearity, the free fermion polarization operator and DSF 
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FIG. 1: Scattering process #2 and bare chiral propaga- 
tors (solid line for right movers, dashed line for left movers, 
wiggly line for the interaction) with non-linear dispersion: 
€t = ±v F(k =F k F ) + (k T k F ) 2 /2m. 




FIG. 2: (a) One- loop polarization diagram and (b) two- 
loop electron self-energy (for a right mover) for the scattering 
process gi- 



(at T = 0) become 

s (0 W = 



rn 
m 



In 



M 2 -(c_) 



21 



(4a) 



LM 2 -(o; + ) 2 J : 

[e(cj - cj_) - b(cj - ,(4b) 



where the two- parametric family appears: w±(q) = 
vp\q\ ± q 2 /2m. The spectral weight now spreads over 
a range of frequencies: cj_(</) < Co> < cj + (g). The latter 
corresponds to a continuum of single (particle-hole) pair 
excitations^ which is responsible for the strong damp- 
ing of the sound mode. From Eq. pb| we see that the 
/-sum rule of Eq. ([3| is now exhausted by the single-pair 
excitations. 

The second approximation made in the TL model is 
that only forward-scattering processes associated with 
small momentum transfers are taken into account. 
The corresponding dimensionless coupling constants are 
noted #4 and g<i for processes involving the scattering of 
electrons of the same and of different chiralities, respec- 
tively, see Fig. [I] for g^. On the other hand, backscatter- 
ing (#i) and umklapp scattering (#3), which are char- 
acterized by large momentum transfers are neglected. 
Under these assumptions, g± and #2 are invariant un- 
der the renormalization group flow and the TL model 
appears as an infrared fixed point of a more general 
class of "g-ology" models. 5 8 It describes a ID metal; 
i.e., excitations are gapless. The difference with the 
3D (Fermi liquid) case may be seerP^ perturbatively 
from the electron-electron scattering rate and quasi- 
particle residue: 7)7 1 = — 9E(^,fc) and z^ 1 = [1 — 
<9!RE(e, h)/de] e= £ k , respectively, where E is the electron 
self-energy; see Fig. [2^b) for the g^ process. The lat- 
ter determine the dressed single particle Green's func- 
tion, G(k 1 e) = l/(e — — E(fc, e)) and the spectral 
function, A(k,e) = — -9G(fc, e). For a Fermi liquid 

(FL), 3E # 2 max{e 2 ,T 2 }/e F and 5RE(e, k) ~ g 2 e, 

implying the existence of well-defined quasi-particles: 
t fl ^ 9 2 max{e 2 , T 2 }/cf, zfl ~ 1 with a Lorentzian 
spectral line-shape. In ID on the other hand, ^E ~ 
-g$msx{\e\,T} and 5RE ~ g\ e ln(max{|e|, T}/e F ). The 
scattering rate is therefore of the order of the energy 
of the particle: ~ g 2 max{|^|, T}, and the quasi- 

particle residue vanishes logarithmically at low energies: 
z tl ~^~92 m ( max {|6c|?^}/ e F), as in a marginal Fermi 
liquid At very weak couplings, #2 <^ 1, fermionic quasi- 
particles survive above an exponentially small energy 



scale: T > ep exp(— l/g^). However, the linear depen- 
dence on energy of the scattering rate implies that the 
spectral line shape is non-Lorentzian with high-frequency 
tails scaling as A(k, e) ~ l/max{|e|, T}, in contrast to the 
Fermi-liquid case. 

Concerning the bosonic correlation function, the com- 
bined assumptions of forward-scattering interactions and 
linear dispersion imply that there is no higher order cor- 
rection to the lowest order diagram of Fig. [2|a). In 
this sense the TL model is exact at one loop. This 
is the essence of the DL (or loop-cancellation) theorem 
which allows for an exact computation of all correla- 
tion functions at arbitrary interaction strength)^ Be- 
yond weak coupling, single-particle excitations disappear 
and the system becomes a non- Fermi liquid. Concerning 
the collective mode, interactions renormalize its velocity, 
v 7^ vp, and bring a prefactor, 7 7^ 1, fixed by the /-sum 
rule in the expression of the DSF, which may then be 
written as Stl(^^q) = j\q\ S(u — v\q\). Bosons therefore 
remain free and the boson-fermion mapping implies that 
there is no relaxation in the TL model. 

As anticipated in the Introduction, relaxation requires 
both dispersion non-linearity and interactions among the 
electrons. The former is irrelevant in the renormalization 
group sense and is not expected to affect the infrared 
properties of the model. But when combined with inter- 
actions among the fermions the DL theorem breaks down 
and multi-loop corrections appear providing a variety of 
collision processes. Presently, no consistent resummation 
scheme exists to take into account these processes in a 
non-perturbative way. In the following we will therefore 
simply focus on the lowest order non-trivial corrections 
beyond the DL theorem and relate such corrections to 
the equilibration of ID interacting fermions. 



III. PERTURBATION THEORY 

We assume that both the interactions among the elec- 
trons and the deviations from linear dispersion are weak 
and take them into account in perturbation theory. For 
simplicity, we focus on spinless fermions, i.e., fully spin- 
polarized electrons. In this case there is no distinc- 
tion between gi and #2 processes. The combinations #4 
and #2 — 9i are then invariant under the renormaliza- 
tion group flow and have to be taken into account. If 
we assume that electrons interact via the realistic 3D 



Coulomb interaction, i.e., V(q) = — e 2 ln(ga), where q 
is the ID momentum transfer and a the width of the 
wire, these coupling constants depend on momentum 
only logarithmically. In the following we will neglect 
this marginal dependence keeping in mind, however, that 
92-91 = [V(0) - V(2k F )}/27rv F ^#2^0, in the log- 
arithmic approximation. Then, to logarithmic accuracy, 
direct and exchange processes involving cancel each 
other and we are left with the momentum-independent 
process #2, see Fig. [I] This considerably reduces the 
number of diagrams that have to be considered, thereby 
simplifying the diagrammatic study while keeping track 
of the most important microscopic processes. 

With these assumptions in hand, the lowest order non- 
trivial diagrams that need to be considered are the three- 
loop diagrams displayed in Fig. [3j For clarity they have 
been separated into three groups with obvious notations: 
"chiral" H±± diagrams as well as a group of "mixed" 
II± T diagrams generally referred to as the Aslamasov- 
Larkin diagrams in the literature. The sum of all these 
contributions yields the lowest non-trivial correction, due 
to both dispersion non-linearity and interactions, to the 



polarization operator: II^ 2 ^ = II Li + II 



t(2) 



2n^_ 2 l, and 

the corresponding DSF. These diagrams have been com- 
puted to the lowest order in 1/ra in Ref. [16]. For clarity, 
we reproduce here the expression of the imaginary part 
of the corresponding retarded polarization operators: 
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where, for simplicity, a numerical coefficient has been 
absorbed in the coupling constant #2 and the thermal 
factor is given by: T — 1/(1 — e 2t^) — 1/(1 — e a/ 9 ). 
From Eqs. ([5| the three-loop correction to the DSF reads 

2 / 9 \ 2 



V F 



(v F q)' 



■T(T;u>,q). (6) 



At zero temperature, J-*(T; uj, q) = Q(uj — v F q) for uj > 0. 
Equation ([6| therefore corresponds to a continuum of 
two-pair excitations which extends beyond the single-pair 
continuum, at uj ^> v F \q\. Thanks to the high-frequency 
cancellation brought by the diagram of Eq. ( 5b ) this two- 



pair continuum has a vanishing contribution to the spec- 
tral weight and the sum rule, Eq. (|3|, is still exhausted 
by the single-pair excitations. This perturbative result 
has limitations as uj — > v F \q\ where non-perturbative ap- 
proaches are requirecP^EO b u t nevertheless agrees with 
similar results obtained by different methods^^. 

On the basis of these three-loop diagrams we can esti- 
mate corrections to the electron-electron scattering rate 
beyond the DL theorem: r _1 = r^ 2 ^ _1 +r^ 4 ^ _1 + ..., where 
r (2)-i = T -i g^T for a thermal excitation as already 
discussed. The next correction is given by the imaginary 
part of the 4-loop fermion self-energy, Fig.[4j Due to con- 
straints imposed by the #2 process the Aslamasov-Larkin 



FIG. 3: Three-loop corrections to the polarization diagram. 



Top: 11^ (and equivalent ly for II ii upon exchanging solid 
and dashed lines). Bottom: 11^. 






FIG. 4: Four- loop corrections to the self-energy of a right- 
moving fermion (and equivalently for a left mover upon ex- 
changing solid and dashed lines). 



diagrams do not contribute to the fermion self-energy. 
The remaining diagrams describe the scattering of two 
(particle-hole) pairs of different chiralities and therefore 
do carry information about the life-time of the excita- 
tions. The latter is estimated to be 
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(7) 



up to a numerical factor. This correction brings back the 
Fermi energy scale but does not follow the usual Fermi- 
liquid-like scaling: T 2 /e F . From dimensional arguments 
we understand the cubic energy-dependence as the only 
possible scaling to compensate for the square of the Fermi 
energy appearing from the polarization part. 



IV. KINETICS AND THERMALIZATION RATE 



Following Refs. |10|12j we now consider the hypothet- 
ical situation where the system is prepared in an out-of- 
equilibrium state characterized by a small temperature 
difference in the distribution of right and left movers. 
The scattering process brings back the system to equi- 
librium. In the following, we give evidence that the cor- 
responding thermalization rate is given by Eq. (J7|. 

A rigorous approach to kinetics should in principle 
be based on the non-equilibrium Green's function tech- 
nique; see, e.g., Refs. [22] for reviews. Here, we shall 
follow the basic steps of the method on a qualitative 
level. This technique allows one to construct a quan- 
tum kinetic equation (QKE) for a generalized distribu- 
tion function: n(e,k,t). For a Fermi liquid the QKE 
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FIG. 5: Comparison between (a) two- loop and (b) one of 
the four-loop diagrams of Fig. [4] Diagram (a) yields a four- 
fermion collision integral. Diagram (b), together with the 
other diagrams of Fig. [4] yields a six-fermion collision integral 
along the cut passing through five particle lines. 

reduces to a Boltzmann equation: d t n k = —I(n k ), where 
n k = n(e = ^ k ,k,t) is the electron distribution func- 
tion and all the relaxation mechanisms are in the col- 
lision integral, I(n). For ID interacting electrons the 
situation is highly non-trivial. The first difficulty is that, 
as the interaction strength increases, quasi-particles be- 
come ill-defined in ID. To circumvent this difficulty, we 
restrict ourselves to the case of very weak interactions 
and not too low temperatures where, as previously dis- 
cussed, quasi-particles still exist. However, the spec- 
tral function is non-Lorentzian: broader on-shell with 
part of the spectral weight transferred to high frequency 
tails. As a consequence one cannot, in principle, neglect 
off-shell contributions: e ^ For thermal carriers, 

we argue that the long time dynamics is dominated by 
the on-shell contribution.^ Focusing on these long time 
scales, which may be the easiest to access experimen- 
tally, we therefore neglect the off-shell contributions. It 
is then well knowrP^that, to lowest order in interactions, 
the collision integral has a four-fermion structure deter- 
mined by the two- loop self-energy diagram of Fig. [5^a) : 
nkn k '- q (l - n k >)(l n k- q ) ~ n k *n k - q (l - n k )(l - n k >- q ), 
where the two terms correspond to "in" and "out" contri- 



butions. Energy and momentum conservation then im- 
pose: £ k +€ k ,_ q = ^k'+^k-q- In 1D tnis yields: q = and 
hence the identical cancellation of "in" and "out" terms. 
We thus recover the fact that (on-shell) two-body colli- 
sions do not contribute to relaxation in ID™ The next 
correction, which is non-zero thanks to dispersion non- 
linearity, is given by the four-loop diagrams of Fig. [4] 
see Fig. [5jb) for a more detailed example. Tedious but 
straightforward calculations then show that the corre- 
sponding collision integral has a six-fermion structure; 
e.g., for the diagram of Fig. [5jb): n k n kl - qi n k2 _ q2 (l - 
n kl )(l-n k2 )(l-n k - qi - q2 )- n kl n k2 n k - qi - q2 (1 -n k )(l - 
n kl - qi )(l - n k2 _ q2 ). The diagrams of Fig. [| yield a 
non-vanishing contribution to the collision integral even 
when the energy-momentum constraint is taken into ac- 
count. On shell, such three-body collisions therefore pro- 
vide the leading contribution to the thermalization rate 
which scales as: g%T 3 /e 2 F , see Eq. 

V. CONCLUSION 



A brief account has been given of the analysis of 
the lowest order non-trivial corrections beyond the DL 
theorem arising from interactions and dispersion non- 
linearity. Focusing on the case of weak interactions 
and spinless fermions, a rigorous derivation of the three- 
loop corrections to the DSF has been reviewed.^ The 
corresponding four-loop fermion self-energy diagrams 
were related to equilibration processes of ID interacting 
fermions. Qualitative arguments allowed us to estimate 
the temperature dependence of the thermalization rate, 
see Eq. | 7[), in agreement with the literature on the sub- 
ject JMnTsee a i so R e f # [24] for recent progress. The dia- 
grammatic approach can be applied to the quantitative 
computation of other equilibri um p roperties of interest. 
In light of recent literature^ ! 12 * ^ the whole approach 
can be extended to, e.g., various interaction potentials, 
to the #4 process™ and, eventually, to the case of spinful 
fermions. 
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From SeSt > 1 and for an energy dispersion Se ~ T, long 
time dynamics correspond to: t ^> 1/T. The off-shell con- 
tribution is dominated by 2-body collisions which affect 
the dynamics on time scales: ~ l/g%T. The on-shell 



contribution corresponds to 3-body collisions which affect 
the dynamics on time scales: ~ e% / 'g^T 3 . The con- 

dition ^> is clearly satisfied in a broad range of 
temperatures: ei?exp(— l/gl) <C T <C £f, where ^2 < 1. 
The #4 process is subject to an infrared catastrophe, see, 
e.g., Ref. [7] for a review, which is cured by curvature cor- 
rections. 



